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Zero modes of massive standard model fermions have been found on electroweak Z-strings. A 
zero mode solution for a massless left-handed neutrino is also known, but was thought to be non- 
normalizable. Here we show that although this mode is not discretely normalizable, it is delta- 
function normalizable and the correct interpretation of this solution is within the framework of the 
continuum spectrum. We also analyze an extension of the standard model including right-handed 
neutrinos in which neutrinos have Dirac masses, arising from a Yukawa coupling to the usual SU(2) 
Higgs doublet, and right-handed Majorana masses. The Majorana mass terms are taken to be 
spatially homogeneous and are presumed to arise from the vacuum expectation value of some field 
acquired in a phase transition well above the electroweak phase transition. The resulting zero energy 
equations have a discrete zero mode. 



I. INTRODUCTION 

The interactions of fermions with solitons has led to 
the discovery of many novel and important field theoretic 
phenomena such as fractional charges and quantum 
anomalies [||-||. Fermionic zero modes on strings |^|^] 
may play an important role in cosmology since they can 
be responsible for string superconductivity |^ and vor- 
tons 1^ . Although the standard electroweak theory does 
not contain any topological defects, "embedded defects" 
such as Z-strings are possible [lO|-p^. Fermionic (quark 
and leptonic) zero modes on electroweak Z-strings have 
been studied in the literature p^-p^ and can provide 
linked configurations of strings with anomalous values 
of baryonic charge and other quantum numbers p^ , p^ . 
Quark and leptonic zero modes also have an important 
effect on the stability of the Z-string ||l7|-[l9|] . 

The work on fermionic zero modes on Z-strings has 
so far considered massless neutrinos. The conclusion of 
these studies has been |5|-[7|, |l^ , |l5| that there are no nor- 
malizable zero energy solutions for massless particles in 
the vortex background. Some attention had been given to 
massive neutrinos before, but in the framework of strings 
in Grand Unified models ||2^,|2^. A reconsideration of 
both massive and massless neutrino zero modes on elec- 
troweak strings is needed. This is the aim of the present 
paper. 

In Section 2, we present the electroweak standard 
model Lagrangian for leptons and find the neutrino equa- 
tions of motion in a Z-string background. We find the so- 
lutions of these equations of motion, and show that a zero 
mode does exist as part of the continuum. In Section 3 
we present the standard extension of the electroweak La- 
grangian with a massive neutrino and derive the neutrino 
equations of motion. We solve the equations of motion 
in two asymptotic regimes (small and large r) and dis- 
cuss the number of well-behaved solutions. We finally 
solve the equations numerically to obtain the zero mode 



solution. 

There are several index theorems in the literature 
which specify the number of zero modes in a given La- 
grangian. The index theorem in |p2| does not apply to a 
massive neutrino case because it assumes that the source 
of all fermionic masses is the interaction with a single 
Higgs field (which is not true for the right-handed neu- 
trino Majorana mass). The generalized index theorem 
in |2^] and |24| also does not apply because it uses an 
ansatz for the spinor components not general enough to 
cover our case. The generalized index theorem in pc| ] 
applies to the massive neutrino with a Dirac mass and 
either a left- or a right-handed Majorana mass but not 
both. 



II. THE STANDARD MODEL NEUTRINO 

We consider first the Glashow-Weinberg-Salam elec- 
troweak standard model, and focus on the interactions of 
a single lepton family. The fermion content is a left- 
handed SU(2) lepton doublet, ^' = (z/i,,e£)"^ and a 
right-handed electron singlet e]^. Electrons acquire their 
(Dirac) mass from Yukawa coupling to the SU(2) doublet 
Higgs field while neutrinos are massless. (There is no 
right-handed neutrino to participate in a Dirac mass, and 
there is assumed to be no left-handed Majorana mass, 
since there is no fundamental Higgs triplet in the model.) 
The Lagrangian of the theory is: 



-X ($t$ - 7y2)2 _ i^^Mi^^^ _ leJ^.i^'D^eR 
+ h' (e^$t^ + ^'^br) . 
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M/;, EE d^w^ - d^wi + g^'^'^w^w: 



(2) 



(3) 



where t" are the Pauh matrices. As is conventional, we 
define 



Af, = sin6l„W^3 



cos 9wB^ 



where 9^ is the weak mixing angle. The fermionic equa- 
tions of motion for Lagrangian are 



(4) 



In the Z-string ansatz, ^ ^ ^ ^ 



(0, (j)) . In cylindrical coordinates (r, 9, z) the gauge and 

= vfir)e'' (5) 



Higgs fields take the form 



0, ee 



where q = a/2 = •y(p+^^/2, and v{r) and f{r) are the 
Z-string profile functions. 

With this ansatz, the neutrino equation of motion is: 



tl^D^i^L = 



(6) 



A representation of the Dirac matrices in cylindrical co- 
ordinates is 



(7) 
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In this representation, a left-handed Dirac fermion has 
the form: = (a, /?, —a, — /3). 

A general solution to equation (^) can be written as a 
superposition of all modes with definite winding number, 



/3 = 



(9) 



Um and (3m are, in general, complex functions of r. We 
are interested in zero energy solutions to equation (||) 
such that all spinor components fall off outside the string 
core (large r) fast enough to be normalizable and are well- 
behaved at the origin (small r). 

After setting kz — 0, we get a set of two recursive 
equations for the coefficients of a and (3: 



-0'm + - 

m — 1 



-Prn = 
m-1 = 



Let us consider the case m = 0. For large r, v ' 
and the closed system of equations is: 

uja^i + I3'n + -/3o = 
r 

—LdPa + a'_i — 



(10) 
(11) 



Eliminating a_i from the system we get Bessel's equation 
for /3: 



1 



1 







/3o +-/?o+/5o(^' 

r r- 

the general solution of which is given by 

/3o(t^?') = AVwJi(wr) + B^/ujYx{ijjr) 



(12) 



(13) 



Here 3\ and Yi are Bessel's functions of the first and sec- 
ond kind, while A and B are constants. We see therefore 
that the fc^ = solutions comprise a continuum labeled 
by Lo. Using the normalization relation for Bessel's func- 
tions: 



1 

Ji{ujr)Jx[Jr)rdr = - cj) . (14) 

LO 



and a similar relation which holds for Yi , we see that the 
solution is properly (5-function normalized. For 7^ and 
r large, ( [l^ ) is just a cylindrical wave /3o ~ e^"^''/^/r. 
To see what happens with the zero mode, i.e. w = 0, 
solution, we take r fixed but large enough for eq. ( [Tl| ) 
to be valid, say r ^ '7"^- Thus, the limit w ^ im- 
plies that LOT — > 0. Using the asymptotic expansion for 
Bessel's function of small argument we get the neutrino 
wave function: 



/3o = 



D 



(15) 



where Z? is a constant independent of and r. 

The neutrino equations for small r and m = are: 



(16) 



1 



= 

r 



where we used v ~ v^r^ and / ^ /or near the origin. 
Eliminating a-\ from the system we get equation for /?: 
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whose solution is given by 



(17) 



(18) 



where C is a constant and Iq is a modified Bessel function 
of the first kind. Thus, (io ~ y^C near the origin, where 
y/uj is introduced in order to have a consistent normal- 
ization. This is in agreement with the result in |1^. 

Po is the only nontrivial solution of the system ([lOD in 
the zero mode limit, i.e. tor — )■ 0. The solutions for all 
other am and (3m are not valid in the zero mode limit 
because of singular behavior at either small or large r 
which makes solutions not even Dirac delta function nor- 
malizable. 

If we look at the solution for Pq as an isolated zero 
mode, it is obviously not normalizable because the nor- 
malization integral diverges logarithmically p^ . But, we 
have seen that this state is actually part of a continuum 
spectrum of the theory and is Dirac delta function nor- 
malizable. Therefore it is a valid zero mode solution. 



The presence of the Majorana mass term, M^, pre- 
vents us from removing the angular dependence from the 
Dirac equations for the neutrino fields using the stan- 
dard procedure Jl3|,[7t. Instead, we must consider a su- 
perposition of all modes with definite winding number. 
Extending the ansatz (O), we write 



(22) 



After setting w = fc^ = 0, we get a set of four recursive 
equations for the coefficients a, /3, 7 and 6 



a = 




-00 ^"li'^)^''''^' 


-iLut-\-im9 


P = 
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-iujt-\-im9 
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-iujt-^-irnO 



7m 7m - iMjijlj^, 
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- h-qfl« 

- h-qJPr, 



(23) 



(to 



m+2 



-6m - iMR5l_.,-^ = - hr]fam^2 



III. BEYOND THE STANDARD MODEL 

Recent evidence from the Super Kamiokande experi- 
ment strongly suggests that neutrinos are very light, 
but not massless. A simple extension to the standard 
model which incorporates light neutrino masses is to add 
to each standard model family an SU{2)l x U{1)y sin- 
glet right-handed neutrino, and to incorporate both a 
Dirac mass for the neutrino and a Majorana mass for the 
right-handed neutrino. The origin of Af^ is a vacuum 
expectation value of some SU (2) singlet Higgs field, ac- 
quired at some energy scale much greater than that of 
the electroweak scale. Mn can therefore be taken to be 
spatially homogeneous. The modified Lagrangian is 



C = C. 



SM 



(19) 



where v'^ = Cv^ and vj^ — (7, (5, 7, 5). The charge conju- 
gation matrix C in the representation we are working in 
is: 



C = i-i^-f = 



'0-10 

10 

-1 

>0 1 



The resulting neutrino equations of motion are: 



(20) 



(21) 



In the Mr = limit, the equations are not recursively 
connected. 

System ( p3| ) decomposes into two subsystems - (/?, 7) 
and (a, S). Let us first analyze the (/3, 7) equations. 
From (^3|) we see that a self consistent set of equations 
consists of four equations for /3m, P-i-m, 7m and 7_i_m. 
The case to = is the most interesting one because we 
already know that in the Mr = case (e.g. electrons), 
zero modes have to = [p^ . Setting to = in the first 
two equations of the system ( p3| ) we get: 

V 

P'o + -Pa = - h'qflo 
r 

PL, + ^^^/3-i = - /177/7-1 (24) 
r 

7o - iMR'jli = - h-qfPo 

7-1 + ;^7-i - iMRj^ = - h-qfP-i 

An analytic solution to the system (|^ could not be 
found, but we can learn something about the structure 
of solutions by looking at their asymptotic behavior. Ne- 
glecting fermion back reaction to the string background, 
V and / are just the Nielsen-Olesen profiles and at large r 
behave as /(r 00) 1, v{r 00) 1. Assuming the 
asymptotic behavior Pq ~ ae"'', 70 ^ be*^ , ~ ce"'', 
7_i ~ de'^ (a, b,c and d are arbitrary complex numbers 
while s, t, u and I are complex numbers, with a nonpos- 
itive real part) we get following conditions: 



s = t = u* = r 



and 



|MR,|V = (|h77p-s2)2 (25) 



The solutions to this fourth order polynomial equation 
for s are given by 
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2\h7^\^ ±\MrW\Mr\^ + A\h7^\^ (26) 



All four of these solutions for s are real. Two of these 
are positive, giving rise to exponentially growing modes 
and two are negative, giving rise to exponentially decay- 
ing modes. Only the latter are physically acceptable. In 
the physical limit of Mr 3> \hri\ we get s+ = Mb_ and 
s- = {hrf)^ / {2M r) ^ which is the standard "see-saw" re- 
sult — outside the string core neutrino zero modes decay 
such that one decay channel is driven by the heavy neu- 
trino (vacuum) mass eigenstate and another one is driven 
by the light one. In the limit of Mr/ {hi]) — > we have 
the standard result s — ±\hr]\ 0. 

If we want to match the exponentially decaying large 
r solutions in a given pair of equations (/3,7) to the solu- 
tions at the origin, we must know how many solutions for 
small r are well-behaved for the set of (/3o, 70, 7-1) . 
By well-behaved we mean nonsingular and single- valued. 
We require not only the functions, but also all their 
derivatives to be single-valued in order to have a well 
defined solution. Therefore, for a solution of the form 
j,Pgiqe ^j^g following must be true: (i) p > \q\, (ii) both 
p and q are even or both are odd. Neglecting fermion 
back reaction on the string background, the asymptotic 
behavior of the profile functions is v ^ v^r'^ + V2r* + . . . 
and / ~ for + f^r^ + ■ ■ ■ for small r. Writing /3o ar^ , 
7o ~ 6r*, f3-i ~ cr", 7_i ~ dr' (a, b,c and d are arbitrary 
complex numbers and s, t, u and I are nonnegative real 
numbers, s and t even, u and I odd) for small r we find 
that the leading orders are s = 0,t = 0,u=l and ^ = 1. 
System ( |2^ ) is invariant under the parity transformation 
(r — > — r), so all corrections to the solutions which are 
parity violating are zero. We can write the general well 
behaved solution in the form: 



1 ^° ^ 




( aor^ - 


h a2r'^ - 


\- a4r* + 




70 








h 64^4 + . 




/3-1 




cqt^ - 


1- C2r^ - 


h C4r^ + . 








\ dor^ ' 


h d2r^ ' 


f diV^ + 





(27) 



The structure of the system ( [2^ ) is such that only three 
coefficients are independent (say Oq, &o and do are inde- 
pendent, while Co = iMRbo/2 and all other coefficient 
are functions of the first three) which says that there are 
three linearly independent well-behaved solutions. For 
the system of four linear first order differential equations 
we should have four solutions in total, so we conclude 
that one solution is not well-behaved, i.e. contains one 
or more terms such as r^, p < and/or r^'(log(r))'. 

With three well-behaved solutions at the origin and 
two well-behaved solutions at infinity there must be 
at least one solution which is well-behaved everywhere. 
Each of the two well-behaved solutions at infinity 
matches on to a unique linear combination of all four 
solutions at the origin where only one of them is bad, so 
there is always one linear combination of the two good 



solutions at infinity which does not have any contribu- 
tion from the single bad solution at the origin. However, 
it could happen that there are two solutions which are 
well-behaved everywhere if each of the two good solu- 
tions at infinity matches exactly to a linear combination 
of the three good solutions at origin. (We call this possi- 
bility a "conspiracy".) One must also be concerned that 
the regular solution for one spinor component, say /3o, 
might correspond to a solution of other spinor compo- 
nents which are not well-behaved. However, it is easy to 
show that if /3o ^ as r ^ 00 then 70, 7-1 — > 0. 
Thus, barring conspiracies, we expect to find only one 
regular solution of the (/?o,7oj/9-i,7-i) equations. 

Numerically solving system (p^, we find a well be- 
haved solution (Figure 1). We took /(r) = tanh(r) and 
v{r) ~ tanh2(r) for the string profile functions which 
gives the correct asymptotical behavior / r, w ~ 
as r — > 0, / ~ 1, w ^ 1 as r cxD and the correct par- 
ity. We also set iMRl(hr\) = 1.1. Although no other 
well behaved numerical solutions were found, a "conspir- 
acy" can not be entirely dismissed because of possible 
numerical errors. However, more than one neutrino zero 
mode would lead to a mismatch between the number of 
neutrino and electron zero modesj^. 




FIG. 1. Numerical solution of the system in eq. ( P4[ ) show- 
ing /3o (solid), 70 (short dash), (dot dash) and 7-1 (long 



dash). The solution is for iMujihr]) 
and v{r) — tanh^(r) 



1.1, /(r) = tanh(r) 



Next we consider non-zero values of the parameter m 
in eq. ([2^). For the case m = 1, analysis of large r gives 
results similar to the m = case — two exponentially 
growing and two exponentially decaying solutions, the 
latter ones with the expected large Mr "see-saw" behav- 



*An analysis in the case of electrons, similar to the one done 
for neutrinos, shows that there is only one well behaved solu- 
tion without the possibility for a conspiracy. 
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ior and the standard Mr limit. For small r, the 
general well behaved solution is of the form: 



1 \ 




( aor^ - 


h a2r^ - 


\- a^r'^ + 


"] 


71 






\- b2r^ - 


h b^r^ + . 








car^ - 


\- C2r** - 


h C4r^ + . 




\l-2 ) 




\ dor^ - 


h dar^ - 


'r d^r^ + 


■■) 



(28) 



where only two coefficients are independent (say bo and 
Co, while ao = -foh-qb^/A, do = iMRbo/A, C2 = -(cqWo + 
hr]foiAI]ibo/4:)/2). With only two good asymptotic solu- 
tions, both at infinity and at the origin, we cannot say 
anything definite about existence of the zero modes in 
this sector. Depending on the specific matching coeffi- 
cients between small and large r solutions we could have 
none, one or two zero modes. However, numerical anal- 
ysis shows that there are no nontrivial well-behaved so- 
lutions, i.e. /3i = 7i = /3-2 — 7-2 — everywhere. 
Numerical experiments also indicates that a similar situ- 
ation occurs for m > 2. 

As a simple consequence of the fact that zero modes 
in 3 -|- 1 dimension satisfy the relationship co = ±/c, they 
must be eigenstates of the 7°7^ operator. We found one 
neutrino zero mode with /3 0, 7 ^ and therefore, 
analysis of the {a,S) pair of the equations of the system 
(^) is not necessary. In order to be eigenstate of 7^7^ 
with the definite eigenvalue (in this case —1) the neutrino 
zero mode must have a = S = 0. 

IV. CONCLUSION 

We have found that in the standard model a massless 
neutrino has one well-behaved zero mode on a Z-string 
which, being part of a continuum, is Dirac delta function 
normalizable. 

We considered the minimal extension of the standard 
model Lagrangian which includes a phenomenologically 
valid massive neutrino, with both Dirac and Majorana 
mass terms. The homogeneous, non-winding Majorana 
mass term A'ln makes this case considerably different 
from the electron one. To accommodate the angular 
dependence of the equations of motion we must use an 
ansatz in which each spinor component is a superposi- 
tion of at least two modes with definite winding num- 
ber. As a consequence, we get a set of coupled equations. 
We solved the equations analytically in two asymptotic 
regimes, small and large r, and argued that in each sector 
determined by a definite winding number m it is possible, 
in principle, to have one or more zero modes. Numeri- 
cal analysis indicates that there is one well-behaved zero 
mode in the m = sector and none in higher m sectors. 
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